Universal asymptotic behavior in flow equations of dissipative systems 
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Based on two dissipative models, universal asymptotic behavior of flow equations for Hamilto- 
nians is found and discussed. Universal asymptotic behavior only depends on fundamental bath 
properties but not on initial system parameters, and the integro-diflerential equations possess an 
universal attractor. The asymptotic flow of the Hamiltonian can be characterized by a non-local 
differential equation which only depends on one parameter - independent of the dissipative system 
or truncation scheme. Since the fixed point Hamiltonian is trivial, the physical information is com- 
pletely transferred to the transformation of the observables. This yields a more stable flow which is 
crucial for the numerical evaluation of correlation functions. Furthermore, the low energy behavior 
of correlation functions is determined analytically. The presented procedure can also be applied if 
relevant perturbations are present as is demonstrated by evaluating dynamical correlation functions 
for sub-Ohmic environments. It can further be generalized to other dissipative systems. 
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I. INTRODUCTION 

The flow equation method, introduced by Glazek and 
Wilsoni and independently by WegnerjS has become a 
powerful tool to analyze many-body systems in nuclear 
and solid-state physics** The method is based on the con- 
tinuous mapping of a given Hamiltonian on to its uni- 
tarily equivalent representations. This induces a flow of 
the system parameters which is governed by the so-called 
flow equations. The continuous transformations can be 
characterized by the anti-Hermitian generator rj and the 
flow equations of the Hamiltonian H are then given by 
the differential form diH = [r),H], £ denoting the flow 
parameter. The objective is to reach a simple fixed point 
Hamiltonian for I — * oo. 

The flow equations for many-body Hamiltonians are 
non-linear, integro-diffcrcntial equations and in general 
analytically non-approachable. Yet, it is possible to an- 
alyze the asymptotic behavior as was done for the spin- 
boson model* as well as for an electron-phonon model. 5 
It was found that the asymptotic behavior of a given 
energy scale A generally leads to A — * A* + a£ -1 / 2 , 
with A* = A(£ = oo) and a denoting a dimensionlcss 
constant. 6 So far, 77 has been chosen such that A* was fi- 
nite and the asymptotic fixed point Hamiltonian has thus 
been characterized by this renormalized energy scale. 

In this paper we want to apply the flow equation tech- 
nique to dissipative systems. The general Hamiltonian 
then consists of a system Hamiltonian H$ with energy 
scale Ao, a bath Hamiltonian Hb, and the interaction 
term Hsb- Instead of integrating out the bath degrees 
of freedom and thus reducing the effective configuration 
space, the flow equation method performs a sequence of 
infinitesimal unitary transformations with the objective 
that the interaction between the system and the bath be- 
comes zero, i.e. Hsb ™* 0. This is done systematically 
by first "rotating away" interaction terms which connect 
states with high energy differences. If A* is finite, the in- 
teraction between the system and the environment is thus 



transformed away in the following way: first decoupling 
the high-energy modes of the bath, then proceeding with 
decoupling the low-energy modes and in the end decou- 
pling nearly degenerate states around the renormalized 
energy scale A r = A* 7^ 0. Whereas in the first two 
regimes the coupling constants decrease exponentially as 
functions of £, the coupling constant belonging to the 
renormalized energy scale A* follows an algebraic decay 
according to £~ 1 ^ 2 , see Ref. 0. 

In order to evaluate spectral properties within the 
flow equation approach, observables need to be subjected 
to the same sequence of unitary transformations. This 
leads to flow equations for observables and the observ- 
able flow follows the systematic decoupling procedure of 
the Hamiltonian, i.e. the correlation function is first built 
up for high and low energies and the intermediate energy 
regime around A* is determined at last. 

Flow equations are more stable for ^-values far away 
from degeneracies. This means that the regime which 
comprises most spectral weight is determined by flow 
equations which are the most unstable. In case of the 
spin-boson model, Kehrein and Mielke circumvent this 
problem by mapping the asymptotic flow of the spin- 
boson model onto the asymptotic flow of the exactly 
solvable dissipative harmonic oscillator £ They thus in- 
tegrated the flow equations numerically up to some finite 
I* and then employed a conservation law which was de- 
rived for the exactly solvable model. Nevertheless, this 
mapping was restricted to a certain parameter regime 8 
and might not be applied to other, more general dissipa- 
tive models. 

In this paper, we will set up flow equations with a 
different asymptotic behavior. The crucial point is that 
the fixed point Hamiltonian H* = H(£ = 00) does not 
contain any characteristic renormalized energy scale, i.e. 
A* = and H* is thus just given by the non-interacting 
bath. As a consequence, the system is monotonously de- 
coupled from the bath, starting from high-energy modes 
and finishing with the low-energy modes. But it also im- 
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plies that the fixed point Hamiltonian does not contain 
any other information but the ground-state energy. If 
the physical renormalized energy scale is denoted by A r , 
the flow equations do not yield A r = A* as before. All 
physical information is transferred on to the observable 
flow. 

The asymptotic behavior of the Hamiltonian is passed 
on to the observable flow as before. Since the asymp- 
totic flow equations deal with the decoupling of the low- 
energy modes, the asymptotic behavior of the observable 
will determine the low-energy regime of correlation func- 
tions. The low-energy regime of dissipative systems is 
generally determined by universal power-laws. We thus 
find the connection between universal flow equations and 
universal spectral properties. We will further set up a 
conservation law, based on asymptotic, scale-invariant 
spectral functions, which will allow us to halt the numer- 
ical integration after some finite value £*. Like this, we 
also obtain analytical results for the low energy behavior 
of correlation functions. The procedure should not be 
restricted to specific models as before. 

There are other advantages of flow equations which 
lead to a trivial fixed point Hamiltonian only consisting 
of the free bath. First, one can easily include a logarith- 
mic grid around uj = to increase the numerical precision 
while decoupling almost degenerate states. Second, the 
spectrum of a dissipative system at low coupling^ is gen- 
erally that of the free bath, consisting of a ground-state 
and a gap-less continuum^fl Unitary transformations do 
not alter the spectrum, but if the fixed point Hamilto- 
nian consisted of the system with the renormalized en- 
ergy scale A* 7^ plus the free bath, the continuous 
spectrum of the bath would be superposed by the discrete 
spectrum of the system. Universal asymptotic behavior 
thus also guarantees the correct spectrum. 

The paper is organized as follows. In Sec. II, we will 
treat the exactly solvable dissipative harmonic oscilla- 
tor. We will first summarize the analytic solution of 
Kehrein and Mielke in a comprised form and then in- 
troduce and discuss universal asymptotic behavior. We 
explicitly show the connection between the conservation 
law obtained from the exact solution and an asymptotic, 
scale-invariant spectral function obtained from the uni- 
versal behavior. In Sec. Ill, we will consider the spin- 
boson model for two different truncation schemes. In 
Sec. IIIA, we investigate flow equations stemming from 
a form-invariant truncation, i.e. no additional coupling 
terms are generated. This will allow a direct comparison 
with the results of the dissipative harmonic oscillator, 
but will yield rather poor results for dynamical correla- 
tion functions. In Sec. IIIB, we will therefore extend the 
truncation scheme and allow an additional coupling term 
to be generated. This approach will yield very good re- 
sults for correlation functions within a large range of the 
parameter space including sub-Ohmic baths. The reader 
who is only interested in this non-trivial application, can 
directly go to Sec. IIIB, which is self-contained. 



II. DISSIPATIVE HARMONIC OSCILLATOR 

We will first discuss the harmonic oscillator coupled to 
an environment. Following the seminal work by Caldeira 
and Leggettjii we will model the bath as a set of non- 
interacting harmonic oscillators with a dense spectrum. 
We will also introduce the interaction induced renor- 
malization of the potential so that the Hamiltonian is 
bounded from below. The dissipative harmonic oscilla- 
tor shall thus be described by the following Hamiltonian: 

^ = y + w+E(f + r«(^-^) 2 ) (i) 

a a 

The parameter k can be decomposed as k — Vo/(77i<7q) 
where Vq shall denote the potential energy scale, m the 
mass of the particle, and go a length scale such that v 
is dimensionlessJ^ We identif y the frequency of the har- 
monic oscillator as ^Jlnv jm. The operators obey the 
canonical commutation relations which read (ft = 1) 

[q,p}=i , [x a ,p a >} = i6 at0l f . (2) 

The system is exactly solvable and has been investi- 
gated by means of several techniques M> The exact solu- 
tion of the model via the flow equation approach was first 
obtained by Kehrein and Mielkeii In Sec. Ill Al we will 
recall the solution in a comprised form in order to intro- 
duce notations. In Sec. IIB, we will then discuss flow 
equations which exhibit universal asymptotic behavior. 
In Sec. Ill CI numerical results are presented. 



A. Analytical Results 

1. Hamiltonian flow 

In order to solve the dissipative harmonic oscillator 
via flow equations, 2 the generator r\ of the infinitesimal 
unitary transformation is chosen to be 

V = i^Q^VlPa +P^fc + E Va,a'X a p a ^ (3) 
a a a, a' 

= if + rf + r\ B ■ 

The first two terms follow from the canonical choice 
rj c = [H, V] with the off-diagonal part V — ~qJ2 a ^c,x a , 
but with generalized parameters rj^, and rf^. They will 
be determined later. The last term rj B is needed to can- 
cel new interaction terms between different bath modes 
which are generated by the flow equations. The solvabil- 
ity of the model is due to the fact that this cancellation 
is possible and that therefore the flow equations close 
exactly. 

The constants of the generator are now chosen such 
that the Hamiltonian remains form-invariant during the 
flow, i. e. no new interaction terms are generated. This 
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procedure is not unique; a possible parametrization is 
given by 

Va = A Q /(w Q ,f) , r\ q a = -X a f((J a ,£) , 

AqAq 



' (f( Wa ,£) + f(u; a ,,£)) 



where f(tv a ,£) denotes an arbitrary function. 

The flow of the system is characterized by the renor- 
malized potential 



2ul 



(4) 



Defining the effective frequency A 2 (£) = 2kv(£) and the 
spectral coupling function 



£ — oo) one imposes 
for £ — ► oo which 



also for any f(ui a ,£). With z* = A 2 
the boundary condition R(z*,£) — > 
guarantees that system and bath are decoupled for £ 
oo. This yields a self-consistent equation for A* = A(£ 
oo) which can be solved for £ — 0. Generally it reads 



A* 



A 2 {£) + A(£) 



J(uj,£)uj 



(10) 



We already want to stress that for the above argu- 
mentation A* has to be finite. Otherwise, Eq. I|l()|) at 
£ = and Eq. cannot hold simultaneously unless 
v(£ = 0) = 0. The case u(£ = 0) = describes the 
system of a dissipative free particle. 

In the next subsection, we will choose f(oJ a ,£) such 
that A* = 0, and the system will still be decoupled from 
the bath for £ — ► oo. The resulting flow equations thus 
belong to a different universality class since Eq. i|l(J|) 
does not hold anymore. This is the key observation of 
this subsection. 



the flow equations d(H = [f],H] are given by the follow- 
ing coupled integro-differential equations: 14 



d e A(£)= / (L,J(oj,£)u>f(oj,£) , (6) 
Jo 

diJ(u,t) = 2J(cu,£)(u J 2 - A 2 (£))f(cu,£) - J( W ,£)M|1 

- 2A(£) J(u, £) [°° du' J{ f' W F(u, J, £) , (7) 
Jo ^ - w 

where we defined F(w,w',£) = f(u,£) + f(w'J). No- 
tice that the renormalization of the bath modes u a was 
neglected since it vanishes in the thermodynamic limit 
which is a typical feature of dissipative systems^ Fur- 
ther, there is no mass renormalization in this approach, 
and we disregarded the renormalization of the ground- 
state energy. 

To solve these equations, Kehrein and Mielke intro- 
duced the function 



R{z,£) = Y J 



A*(*) 



z — w± 



(8) 



for which the following differential equation holds: 

diR(z,£)^~diA 2 (£) (9) 
+ 2(z- A 2 {£) - R(z, £)) ^i/k^) 

a u 

The algebraic equation z — A 2 (£) — R(z,£) — solves 
the above differential equation for any z, and apparently 



2. Observable flow 

To determine correlation functions we have to apply 
the same sequence of infinitesimal transformations to the 
observables that led to the diagonalization of the Hamil- 
tonian. Only then the diagonal structure of the Hamilto- 
nian can be used to yield a simple time evolution of the 
operators in the Heisenberg picture. For the flow of the 
position operator we make the ansatz 



q(£) = h(£)q + J2xa(£)x a 



(11) 



The initial conditions are given by h{£ = 0) = 1 and 
Xa(£ = 0) = 0. During the flow, the weight of the system 
operator will be transferred to the bath operators. In 
the language of the flow equation approach we speak of a 
dissipative system when the total weight of the system is 
being transferred to the bath during the flow, i.e. h(£ = 
oo) = 0. The flow equations for the observable diq = 
[rj, q] close and with the spectral function 



S(cj,£) = -J2 



K(£)Xa(£) 



S(lu — U! a ) 



(12) 



we obtain the following integro-differential equations: 

d e h(£)= duL>S(u,£)f(u>,£) , (13) 
Jo 
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d e S(oj,£) = -A(£)h(£)J(oj,£)f(oj,£) - A{£)J{oj,£) 
+ (oj 2 - A 2 (£))S(oj,£)f(oj,£) - A(£)S(oj,£) 

I 



dw — ^-F(w,w ,£) 

o oj- — OJ 



The flow equations for p(£) = h(£)p + J2 a Xa{£)Pa are 
equivalent to Eqs. (|13l) and (|14H which merely demon- 
strates the fact that Hermite-city is conserved during the 
flow. From the above flow equations we obtain the fol- 
lowing sum rule, which signifies that the commutation 
relation [q(£),p(£)] — i holds for all £: 



(15) 



To solve the flow equations one introduces the functions 

2 ' 

a 

(16) 



S 2 (z,£) = J2 



/(a 



Kehrein and Mielke showed that the following quantity 
is conserved: 



-&(*,£)) 2 
S ^ t)+ z -A(£)-R(z,£) = C ° nSL 



(17) 



doj — ^-F(oj,oj ,£) 

0J Z — oj' 



(14) 



One can show that h(£ = oo) = for the initial conditions 
of interest such that the total weight of the system is 
indeed transferred to the bath. Correlation functions are 
obtained through the following identity: 



K(u;,£)^J2^ 6 ^ 2 -^ 



= -lmS 2 {oJ 2 - i0,£) 



(18) 



For example, the spectral function ojK(lo) = ojK(oj, £ 
oo) is proportional to the Fourier transform of (q(t)q). 

Choosing a Lorentzian spectral function 

Arfoja. 



J(u) = J(u,t = 0) „ , 
or + 7 Z 



the explicit solution is given by (uiq = A(£ — 0)) 7 



(19) 



K(uj) 



2a7 2 w w(7 2 + uj 2 ) 



(oj 2 (j 2 + oj 2 ) - 2iraj 3 oj Q - oj 2 (-f 2 + oj 2 )) 2 + 4ir 2 a 2 oj 2 j 4 oj 2 

I 



(20) 



B. Universal Asymptotic Behavior 



In the previous subsection, the generator of the flow 
equations, rj, was chosen such that the initial Hamilto- 
nian of Eq. Q remained form-invariant. This condition 
defines the generator only up to an arbitrary function 
f(oj,£). It is often useful to explicitly specify f(oj,£) 
even though the final result must be independent of 
the particular choice. In Ref. [jj, Kehrein and Mielke 
choose f(oj,£) = -(oj-A(£))/(oj + A(£)). This leads to 
J(oj, I) ex exp[— 2(oj — A(£)) 2 £] for £ — > oo if one neglects 
the non-linear term in Eq. l[7)l. The spectral function 
is therefore centered around the renormalized frequency 
A* 7^ 0. At uj = A* the spectral function vanishes al- 
gebraically as £~ 1 / 2 . The asymptotic spectral function 
thus depends on the initial frequency through A*. We 
want to label this asymptotic behavior non-universal. 



A different choice is f{oj,£) = —1 which would have 
the consequence that the renormalized potential has to 
tend to zero so that the system is decoupled for £ — > 
oo, i.e. A* = in order that J(w,£ = oo) = 0. It 
therefore belongs to a different universality class since Eq. 
HI U| l does not hold anymore. The differential equation 
@ turns into a Ricatti equation which can formally be 
integrated. 

That the choice f(oj,£) = —1 really decouples the sys- 
tem from the bath is not clear from the beginning. One 
way to convince oneself is to study the asymptotic be- 
havior of the flow equations ijfjj and Q. For this, we 
introduce an additional energy scale related to the cou- 
pling function, A(£) = / °° dojJ(oj, £)/oj. 

We now make the ansatz 



A(£) -> al 



-1/2 



A(£) -> b£ 



-1/2 



(21) 



with constants a and b as 



oo. We further assume 



5 



tha 
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j(cj,e) -» lo s ' j{i)j(uj^Ti) 



(22) 



The various contributions of J(w, ^) can be interpreted 
as follows: The energy dependence uj s will yield the (cor- 
rect) algebraic low-energy behavior of the spectral func- 
tion, J(ujy£) represents the high-energy cutoff function 
with the ^-dependent cutoff frequency lo 1 c — £~ x l 2 , and 
the function J(l) is needed to assure units of energy. 

The differential equations for K(£) and J(£) can be 
obtained from Eq. J7J) and the special feature of the 
exactly solvable dissipative harmonic oscillator is that the 
differential equation for A(£) closes, i.e. only A(£), A(£), 
and diA{£) are involved. One obtains J(£) — > £^ s _1 )/ 2 
and 



= l/2 + s'/4 , 6 = a/2 



(23) 



We now define the dimcnsionlcss, scale- invariant vari- 
able y = ujVI and set J(y) — > Jo for y — > 0. We then 
obtain the following non-local differential equation for 
J(y) ee V s '- 1 J{y): 



dyJ(y) 



4yJ(j/)(l- 
J(y) 



2 a 



dy'- 



Ay') 

yl _ yl- 



(24) 



The boundary condition is given by a = 2 dyJ{y) . 

The above ansatz for the asymptotic behavior (Eq. 
(|2f|l and H22[)) guarantees that the system will be de- 
coupled from the bath since the support of the spectral 
function vanishes as t~ x l 2 and J(y) — > y 4+ ^ s ~^e~ 2y ~ for 
y — > co. There is thus a universal fixed point for all initial 
frequencies. They are all mapped onto the free particle 
plus bath. 

To determine the asymptotic behavior of h{€) and 
S{u),£), we need the quantity T,(£) = J °° duiS(u},£)/uj. 
We make a similar ansatz as in the case of the spectral 
function, namely 

h{£) -> cT 1/2 ^ , ^{t) ^ d£- 1/2 -t (25) 



with constants c and d as 



and 



S(w,*)->w* S(£)S(y) 



(26) 



Notice that we introduced a further parameter £. 

The differential equation for E(£) and 5(f) are ob- 
tained from (HU). One obtains S(£) -> £-Z+{s'-i)/2 an d 



= 2d 



(27) 



Setting S(y) — * So for y — > 0, we obtain the following 
non-local differential equation for = y s _1 5(y): 



5,5(2/) = -2yS(y)(l -2a I" dy'^l^j 



2 , (28) 



+ 4ayJ(y) [°° dy'-p^ + (s' - 1) 

■/o ?/ 2 - y' 2/ 



Eq. I|28[) is linear in S(y). We can therefore not determine 
the constant So from the asymptotic behavior. But we 
can extract the low-frequency behavior of the correlation 
function K(ui) = lim^oo S 2 (u;, t)/{2A{£) J(u, £)): 



K{lo -^0)=lo s 



o2 
°0 

2aJn 



(29) 



The numerical results yield s' = s, where s denotes the 
algebraic behavior of the initial spectral function at low 
energies, i.e. we recover the well-known result. 

Finally, we will give an alternative derivation of the 
conservation law of Eq. I|f 7[) which is valid in the asymp- 
totic regime. For this, we define the asymptotic function 
K a (u>,£), following from the conservation law: 

K a {u,£) = — — — — — (30) 

7r uj z — A z (£) — R(Lu z — iv, £) 

Taking £* to be in the asymptotic regime, we obtain 
K a (u,£) ee K(uj) - K{u,t) = uj s ' ' K a (y*) with y* = 
uj^/I* and 



K„ ( a . = 4 / dyy 2 - s 'S{y) [°° dy' 

v* Jo y y 



(31) 



One can thus derive the conservation law from the uni- 
versal asymptotic behavior, given that £* is in the asymp- 
totic regime. This is an important observation since con- 
servation laws as Eq. I|f 7fl can only be derived for solvable 
models. But by analyzing the asymptotic behavior, we 
are able to integrate the flow equations up to a finite but 
asymptotic £* and than complete the spectral function 
by analyzing the non-local differential equations for the 
scale-invariant spectral functions. This procedure also 
works for non-trivial models. 



C. Numerical Results 

We will now analyze the above flow equations numeri- 
cally. Given the Lorentzian of Eq. I|f 9fl as initial coupling 
function, we will calculate K(lo) within the flow equation 
approach for two different choices of f(u>,£). To do so, 
we will employ the conservation law of Eq. Q 7fl . i.e. 
K(u) = k\uj,£) + K a (u,£). This allows us to halt the 
integration routine after a finite i*. Nevertheless, the 
final result must be independent of £* which is the case. 

In Fig. ^ the results are shown for a = 0.1 and 
"f/uJo — 1- The analytic solution superposes the 
solution obtained by flow equation approach and, 
obviously, K(uo) = K(cu,£) + K a (uj,£) is indepen- 
dent of the two choices of f(ui,£). Nevertheless, the 
functions K(oj,£) and K a (ui,£) are different. For 
f(u,£) = -(u - A{£))/{u + A(£)), K a (cj,£) mostly 
contains the weight around intermediate time scales, 
whereas for f(uj,£) = —1, K a (ui,£) determines the 
low-frequency behavior. 
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to/ro (o/w 

FIG. 1: The spectral function K(u) = K(u,£*) + K a (uj,t) 
calculated via flow equations with f(u,£) — — (uj — A(£))/(u> + 
A(£)) (lhs) and with f{u,l) = -1 (rhs), taken at Jot* = 10 
for J(uj) = 47 2 wa/(w 2 + 7 2 ) with lo = A(£ = 0), a = 0.1 and 
7/ljo = 1 (dotted line). The solid line resembles the analytic 
solution of Eq. 12011 which is almost entirely superposing the 
dotted line. 



0.5 




FIG. 2: The asymptotic functions J(y) and S(y) at vol* — 
100 for Ohmic coupling J{u) = 2awO(o» c — w) with lj 2 /vo = 
100 and no = kv(£ — 0) for different coupling strengths a. 
The solid lines resemble the "analytic" solution. 



We will now investigate the asymptotic behavior of the 
flow equations for f(u),£) = —1. The differential equa- 
tions (SJl and (j2EJ can be solved numerically via self- 
consistent iteration. Comparing these "analytic" results 
with the numerical results obtained via integrating the 
flow equations shall serve as a proof that we have found 
the correct asymptotic behavior. 

In the following, the initial spectral coupling function 
shall be given by J(u>) — 2aK 1 ~ s ui s O(ui c — iS) where a 
denotes the coupling strength and lo c the high-frequency 




y 



FIG. 3: The asymptotic functions J(y) and S(y) with s' = 1 
at v £* = 100 for J(w) = 2aK 1 ~ s lj s Q{lo c - u) with a = 0.4, 
K 2 /vo — 1, lo 2 /vq = 100 and no = kv(£ = 0) for different 
bath types s (s' — s). The solid lines resemble the "analytic" 
solution. 




FIG. 4: The universal asymptotic function K a (y)/S 2 for sub- 
Ohmic (lhs), Ohmic, and super-Ohmic baths (rhs) (s' = s). 



cutoff. The bath type is labeled by s and K denotes an 
additional energy scale present for s / 1. One distin- 
guishes the following three characteristic coupling types: 
super-Ohmic coupling (s > 1), Ohmic coupling (s = 1), 
and sub-Ohmic coupling (s < 1). The numerical data 
shows that s' = s. 

The results for Ohmic coupling s = 1 are shown in Fig. 
13 where the asymptotic curves J(y) and S(y) obtained 
through direct integration and through self-consistent it- 
eration are plotted for different coupling strengths a. 
Indeed, J(y) is universal for all coupling strengths, the 
functions S(y) only differ by the constant So. This be- 
havior also holds for s ^ 1. 

In Fig. |3 the scale-invariant functions J(y) and S(y) 
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are shown at fixed coupling strength a — 0.4 for two 
different parameters s. In the super-Ohmic case s = 1 .5, 
J{y) and S(y) show a suppression whereas in the sub- 
Ohmic case s — 0.5, J(y) and S(y) diverge as y — > 0. 

We close this section with a discussion on the asymp- 
totic function K a (y) of of Eq. (|31fl which is needed for 
the alternative conservation law. K a (y) is proportional 
to Sq and thus depends on the coupling strength a. Scal- 
ing out this factor, we obtain universal functions only de- 
pending on the bath type (s' = s). In Fig. the results 
are shown for sub-Ohmic (lhs), Ohmic, and super-Ohmic 
baths (rhs). Notice that only for Ohmic dissipation s = 1, 
d y K{y) for y -> 0. 



III. SPIN-BOSON MODEL 

The most prominent dissipative quantum system is 
the spin-boson model. It has been extensively stud- 
ied over the past thirty years by several varying tech- 
niques - ranging from the Feynman- Vernon influence- 
functional formulation^ Liouville operator and projec- 
tion methods^ to exact mathematical results^ Its pop- 
ularity stems from the fact that it is the most simple 
non-trivial dissipative quantum system, consisting only 
of a two-state system and a bosonic bath that are lin- 
early coupled in order to preserve the reflection symmetry 
of the system. Nevertheless, applications for the model 
are found in all fields of physics, starting from quan- 
tum optics and solid state physics to nuclear physics and 
chemistry 13,19 In addition, from a mathematical and the- 
oretical point of view, it attracts a lot of interest since it 
exhibits a Kosterlitz-Thouless transition for Ohmic cou- 
pling at T = and a crossover from coherent to inco- 
herent tunneling at finite temperature for all coupling 
types. Moreover, it is related to other prominent models 
of theoretical physics - most strikingly to the anisotropic 
Kondo model. 

The Kondo model as well as the spin-boson model 
belong to the so-called strong-coupling problems which 
make the use of renormalization group techniques almost 
indispensable. The low-energy behavior of the Kondo 
model was thus first obtained by Wilson employing his 
numerical renormalization grouptSS These ideas were also 
applied to the spin-boson model^i Flow equations results 
for the s pin-b oson model and Kondo model can be found 
in Refs. I7l22ll2l2l 

The Hamiltonian of the symmetric spin-boson model 
without bias is given by 

^ = -y^ + E^ + ^Ey^ + 6 «)- (32) 

a a 

The operators ba resemble the bath degrees of free- 
dom and <7i with i = x,y,z denote the Pauli spin ma- 
trices. They obey the canonical commutation relations 
and the spin-1/2 algebra, respectively. The coupling con- 



stants X a only enter via the spectral coupling function 
JH = V](A»)>-, tt ), (33) 

a 

which shall be parametrized as 

J(lo) = 2aK 1 - s Lu s Q(Lu c - u). (34) 

The dimcnsionlcss parameter a denotes the coupling 
strength and lo c the high-frequency cutoff. The bath type 
is labeled by s and K denotes an additional energy scale 
present for s =/= 1. 

The Hamiltonian l)32|l resembles an effective Hamilto- 
nian where the high energy degrees of freedom of the 
bath were already integrated out - down to the bath cut- 
off u> c by employing the Born-Oppenheimer approxima- 
tion. The tunnel-matrix element Ao thus depends on 
this arbitrary energy scale and all physical results must 
be independent of the cutoff u> c . For Ohmic coupling one 
obtains A cx (w c ) a and the only combination of A and 
lu c that yields the units of an energy and is independent of 
uj c is given by A r cx A (Ao/uj c ) a ^ 1 ^ a K Physical observ- 
ables should only depend on this effective tunnel-matrix 
element. 

In the following, we will discuss the symmetrized equi- 
librium correlation function 

C(t)= 2<{*.(*W> = y dte wt C(uj) (35) 

at low coupling. In this regime, C(u>) has been eval- 
uated for super-Ohmic baths and for Ohmic baths by 
Kehrein and Mielke employing flow equations which ex- 
hibited non- universal asymptotic behavior Here, we will 
set up and investigate flow equations which exhibit uni- 
versal asymptotic behavior. 

In Sec. IIII Al we first discuss flow equations which pre- 
serve the form of the Hamiltonian of Eq. i|32|) . A direct 
comparison with the flow of the dissipative harmonic os- 
cillator is then possible, but this truncation scheme will 
yield rather poor results for C(u). In Sec. IIII Bl we will 
therefore extend the truncation scheme and also allow ad- 
ditional coupling terms to be generated. This procedure 
will yield very good results for equilibrium correlation 
functions even for sub-Ohmic baths. 

This extended scheme includes operators which we 
want to label "transitionally relevant" . Those operators 
neither appear in the initial nor in the fixed point Hamil- 
tonian, but are obviously needed as transition link. This 
extents the classification of operators for the flow equa- 
tion approach, initiated by Kehrein and Mielke in Rcf. 
l25l We want to note that this "transitional relevance" 
can already be seen from a simpler model including only 
one bosonic modern^ 
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A. Form-invariant flow 

1. Flow equations 

Applying flow equations to the spin-boson model we 
will first be guided by the premise that the flow shall not 
generate new interaction terms. This cannot be done ex- 
actly as in the case of the dissipative harmonic oscillator 
but one has to neglect normal ordered coupling terms of 
higher order. Then, one can follow the same procedure 
as outlined in the previous section. We therefore just cite 
the resul ting equations. For details, we refer to Ref. 7 
and Ref. 

The flow equations diH — [i], H] for the spin-boson 
model with form-invariant flow are given by the following 
coupled integro-differential equations (T = 0): 

/>OC 

d e A(£) = A{£) duJ(u,£)f(u,£) (36) 
Jo 

d e J(cj, £) = 2J(lj, £){u 2 - A 2 {£))f{u J , £) (37) 
- 2A(£) J(u>,£) rdw , J{ f J) ^ F(u,u',l) 

JO UJ Z - LJ' 

We introduced the ^-dependent spectral coupling func- 



where we defined S(uj,£) = J2 a X a (£)xa(£)S(uj — w a ). 

In Ref. |7J, the flow equations are numerically inte- 
grated, choosing f(u),£) = -(u-A(£))/(l> + A(£)). The 
asymptotic spectral function is then centered around the 
renormalized tunnel-matrix element A r and one is able 
to map the asymptotic flow equations of the spin-boson 
model on to the asymptotic flow equations of the dissipa- 
tive harmonic oscillator. One can then employ the exact 
solution of the dissipative harmonic oscillator outlined in 
Sec. IIA. The result is shown in Fig. |5J 



2. Asymptotic behavior 

In analogy to the dissipative harmonic oscillator we 
now investigate the flow equations with f(ui,£) = — 1 
which will lead to universal asymptotic behavior. These 
flow equations show apparent deficiencies when calculat- 
ing correlation functions. E.g., the sum rule (<J 2 (£)) = 1 
is not satisfied for f(ui,£) — —1 which can be seen from 



tion J(uj 7 £) = Ea^aW^( w — w a) an d defined 
F(w,uj',£) — f(to,£) + f(uj',£). Further, we disregarded 
the renormalization of the ground-state energy. 

As in the previous section, the condition of form- 
invariance does not specify the flow completely and the 
flow equations still depend on an arbitrary function 
f(u>,£). But if one considers the spin-boson model with 
finite bias or starts from a unitarily equivalent repre- 
sentation of the initial Hamiltonian where trje reflec- 
tion symmetry is broken, e.g. b a — > b a — ^-f-, there 
is no open choice for f(ui,£) anymore and f(u),£) = — 1 
emerges automatically, 27 Breaking the reflection symme- 
try of the system thus also breaks the "gauge- invariance" , 
expressed by the function f(u>,£). 

The truncated flow of the observable shall be given by 

a z (£)=h(£)a z + a x J2xa(i)(b a + bi) ■ (38) 

a 

The flow equations for the observable di<j z — [77, a z ] are 
approximated to yield 

poo 

d t h(£) = A{£) / dujS(uj,£)f(uj,£) , (39) 
Jo 



(40) 



I 

the differential equation of (al(£)} 7 

dt{ai(£)} = -2A(£) r duriJ 8 ^ 8 ^^ . (41) 

Jo UJ + LU' 

Since h(£) — ► for £ — > oo and S(u>,£) > 0, the spectral 
function C{u>) = J2 a Xa(^ = oo)S(u> — u> a ) cannot satisfy 
the spectral sum rule J °° dcuC(uj) = 1. 

Analyzing the asymptotic behavior, we now show that 
the flow equations do not yield the right asymptotic be- 
havior, either. We limit ourselves to Ohmic dissipation. 

Following the discussion of the previous section, we 
make the ansatz 

A^^ar 1 ' 2 , A(£)^b£-^ 2 (42) 

with A(£) = J °° dui,J(uj,£)/uj as £ — > oo. We further 
assume that J(cj, £) — ► uJ(y), where J(y) — > Jo for y — > 
and y = uVI- We observe that a — 2b and thus obtain 
the following non-local differential equation: 

d y J{y) = -AyJ(y) (l - 2a f°° dy' '-^s) (43) 

v Jo y-y'> 



i 

d t S(u, £) = -A{£)h(£)J(u, £)f(w, £) - A(£) J(«, £) f duj / S ^ F(u, u', £) 

Jo u - oj' 

+ (oj 2 - A 2 (£))S(oj,£)f(uj,£) - A(£)S(oj,£) f°° dJ J(fa/ ' ^ F(oj, oj', £) 

JO ^ - U3 



The boundary condition is given by 1 = 2 J °° dyyJ(y). 

The above equation is identical to Eq. I|24ll with the 
only difference that the parameter a = 2 J °° dyJ(y) takes 
on a different value and cannot be determined analyti- 
cally, now. In the next subsection, we will see that the 
above equation also describes the asymptotic flow of the 
coupling function for an improved truncation scheme, but 
with a = 2. 

To determine the asymptotic behavior of S(lu, £) we 
make the ansatz 

h(£) -> c£- 1/2 - £ , E(e)-+ d£- 1/2 -* (44) 

with T,{£) = J °° dioS(u),t)/uj as I — > oo. Further we 
assume that S(uj,£) — > uS(£)S(y), where S(y) — > So for 
y — > 0. This leads to the following non-local differential 
equation for S(y): 



dyS(y) = -2yS(y)(l-2a 
+4ayJ(y) H ' dy' ' S( ' /) 



dy' 

y^~y 

S(y) - J(y)S /J 



(45) 



2£- 



There is one difference to the corresponding differential 
equation of the dissipative harmonic oscillator since the 
last term in Eq. I|45l) is not present in Eq. I|28|l . 

To see the deficiency of the above asymptotic behavior, 
we now observe that 



L0t 



02 
-2£ J 

Jo 



(46) 



We would only obtain the correct asymptotic behavior 
C(ui) oc lo for u> — ► if £ = 0. As we will see from the 
numerical results, this is not the case. We thus obtain 
the wrong asymptotic behavior which is an evidence for 
the general shortcoming of the form-invariant flow. 

We want to close with a remark on non-universal 
asymptotic behavior. For Ohmic coupling it yields 
A(£) — > A* + \£~ 1 I 2 for coupling constants a < 1, 
with A* denoting the renormalized tunnel-matrix ele- 
ment A* tx A (A /w c ) Q /( 1 - Q ), see Ref. H For a > 1 
the asymptotic behavior is governed by A — > £~ 1 / 2 /ln£. 
The localization phenomena can thus be detected by the 
different asymptotic behavior of A(^). This is in con- 
trast to the universal asymptotic behavior where for all 
coupling constants A(^) — > al^ 1 !" 1 holds. This indicates 
that our approach will only be valid for small coupling 
strength, away from the broken symmetry phase. 



3. Numerical results 

The predictions of the above ansatz can be verified nu- 
merically by integrating the flow equations. In contrast 
to the exactly solvable dissipative harmonic oscillator, 
the constants a and £ cannot be determined analytically. 






Oa = 0.5 




n □ a = 0.4 _ 




\ Ca=0.3 


o * J 


1 







FIG. 5: The asymptotic spectral functions J(y) and S(y) 
obtained from the form-invariant flow for Ohmic coupling at 
Aq£* = 10 3 for different coupling constants a. The solid line 
resembles the "analytic" solution. 




FIG. 6: The spectral functions C(u>) obtained from the form- 
invariant flow for Ohmic coupling J(u>) — 2aujQ(ui c — uo) 
with u! c /Ao = 30 at Aq£* — 10 3 for different coupling con- 
stants q. Inset: Defining the maximum of C(lo) as A r , 
ln(A r /Ao) is plotted versus x(a, ui c ) = a/(l — a) ln(a» c /Ao) — 
ln(cos(7ra)r(l - 2a))/ (2 - 2a) for different cutoffs uo c . The 
solid line resembles the NIBA result for A e ff. 



The self-consistent solution of the differential equations 
iffify and (|^5|) yields the following values: 



1.21 



i w 0.41 



(47) 



In Fig. the universal asymptotic functions J(y) 
and S(y) obtained from numerical integration of the flow 
equations are compared with the "analytic" solution fol- 
lowing from the self-consistent evaluation of Eqs. I|43|) 
and l|45|) . The asymptotic results do only depend on 
the bath type (not shown here), but not on the coupling 
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strength a, the explicit form of the cutoff function, or on 
the value of ui c . 

The symmetrized equilibrium correlation function 
C(w) must of course depend on the initial conditions. 
Since the asymptotic functions S(y) only differ by at 
most a constant for arbitrary initial conditions, this in- 
formation must be contained in the primary and inter- 
mediate Grange of the flow. In Fig. [5] the flow equation 
results of C(uj) for a Ohmic bath with sharp cutoff at 
u> c / A = 30 at different coupling strengths a are shown. 

The spectral function C(u>) is obtained by integrating 
the flow equations up to Aq£* — 10 3 where one is already 
in the asymptotic regime for the coupling strengths con- 
sidered. For higher coupling strengths one would have 
to integrate up to larger £* , since A r is smaller. Notice 
that we do not employ a conservation law to build up the 
spectral functions, nor is our approach limited to small 
coupling strengths. But one crucial shortcoming of the 
flow equations is that C(cu) does not fulfill the sum rule 
J °° du)C(uj) = 1. Independent from the cutoff parameter, 
u c and bath type the sum rule only yields approximately 
80% for small coupling constants (a<0.1) and 83 — 85% 
for larger coupling constants. 

Nevertheless, C(u>) contains physical information. 
This is demonstrated in the inset of Fig. where the 
maximum of C(u>) is plotted for different cutoff frequen- 
cies uj c as a function of the coupling strength a, according 
to the effective tunnel-matrix element 

A cff = ( cos(7ra)r(l - 2a)) ~ 2(1 - q) A (^) ~ , (48) 

derived from the non-interacting blib approximation 
(NIBA) . 19 The universal regime described by the NIB A is 
only reached for small coupling constants and for rather 
large oj c . Still, the flow equations capture the essential 
physics despite of the apparent shortcomings of the cor- 
relation function. 



B. Extended flow equations 

1. Flow equations 

Since the flow equations of the form-invariant flow did 
not yield satisfactory results, we will now extend the 
truncation scheme, i.e. we will choose a different gen- 
erator and we will allow a new coupling term to be gen- 
erated. But universal asymptotic behavior shall be re- 
tained. 



The anti-Hermitian generator rj shall first be chosen 
canonically, i.e. r/ c — [Hb,H], with the diagonal Hamil- 
tonian Hb — J2 a w a^a and H given in Eq. 1)32(1 . This 
choice indicates that there will be universal asymptotic 
behavior since the diagonal Hamiltonian consists of the 
bath only and is thus ^-independent. 28 

The commutator [r) c , H] will now generate a new cou- 
pling term which is linear in the bosonic operators. Ne- 
glecting the shift in the ground-state energy, we thus ar- 
rive at the following truncated Hamiltonian: 

H(£) = -^a x + J2"<*bib a (49) 

a 

a. a 

with A{£ = 0) = A , X z a {£ = 0) = A°, and WJl = 0) = 0. 

The canonical generator r\ c also gives rise to coupling 
terms which are bilinear in the bosonic operators. These 
terms will not be included in the Hamiltonian flow ex- 
plicitly. Still, they are adequately taken into account of 
by introducing an additional term to the generator which 
approximately cancels these bilinear contributions. This 
is in analogy to the procedure of the previous section. 
For details, see Ref. [2J. 

In terms of the spectral coupling functions J k (ui,£) = 
Ea^a) 2 ^)^ — wrtn k = y,Z, we obtain the fol- 
lowing flow equations: 

/>oo 

diA = -2 dujujJ m (uj,£) (50) 
Jo 

deJ k (uj,£) = ~2uj 2 f k (uj,£) + 2AujJ m (uj,£) (51) 

poo 2 I /2 

+ 2J k (w,£) / dw'J m (cj',£) ^ U /2 
Jo u - v 

+ U m (cj,£) / do/JV,*) r ,2 
Jo w - 

We defined k = (z,y) for k = (y,z) and J m (uj 7 £) = 

J2 a KK*(" - "a) = (J»(w,*)J*(w,*)) 1/2 . 

The truncated flow of the observable shall be given 
by a z {£) = h{£)a z + a x J2 aXa (£)(b a + ftt). The flow 
equations di a z — [77, a z ] are subjected to the same de- 
coupling procedure as above. Defining the spectral func- 
tions S k {uj,£) =E a ^oWXaW^(w-w a ), we obtain the 
following flow equations: 

poo 

deh(£) = - (L>uS v (luJ) , (52) 
Jo 
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d e S y ((J, I) = h(£)ujj y (uj, £) - uj 2 S v (uj, I) + A(£)ujS z (uj, £) 

r°° ll> 2 + ui' 2 r°° 

+ J m {u,£) / dui'S y (uj',£) 2 + 2J y (uj,£) / duj'S z (uj',£) 
Jo oj 2 - oj' Jo 

r°° uj 2 4- uj' 2 r°° 

+ S z {ujJ) / duj'J y (uj',£) — - + 2S y (u>,l) / duj'J m (uj',£) 
Jo uj 2 - uj' Jo 



uj 2 - uj' 2 



uj 2 - u' 2 



(53) 



Since S z (oj,£) = S y (uj,l)J z {uj,£)/J m (uj,l), it suffices to 
define the integro-differential equation for S v (uj,£). 

One criterion for the assessment of the quality of the 
flow equations is again given by (<J 2 {£)) which should be 
approximately equal to one for all I. It is governed by 
the differential equation 

/•oo / 

dAa 2 J£))=- dujduj'— — S y (uj,£)S z (uj',£). (54) 
Jo uj + uj' 

In the asymptotic regime I — ► oo, the flow equations 
simplify considerably, i.e. S v (uj,£) — S z (uj,£) which 
will be shown below. This property guarantees that the 
spectral weight of correlation functions is asymptotically 
conserved - in contrary to the form-invariant flow. Nu- 
merical calculations show that - throughout the flow - 
the spectral sum rule is fulfilled within less than 0.1% 
relative error. This is a significant improvement w.r.t. 
the form-invariant flow. 



2. Asymptotic behavior 
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FIG. 7: Left hand side: The asymptotic spectral functions 
J(y) obtained from the extended flow for Ohmic coupling with 
a — 0.1 and lj c /Ao = 10 at various values of i. Right hand 
side: The "analytic" solution following from Eq. 1571 for the 
corresponding values of a(l) = A(^)^ 1//2 as they follow from 
the flow equations. A numerical fit suggests a — > 4.5 for 
t -> oo. 



We now discuss the asymptotic behavior. As a first 
approximation, we take A* = \ z a and thus drop the index 
k of the spectral function. This assumption is consistent 
with the flow equations and justified later based on the 
numerical result. The asymptotic flow equations are then 
solved by the ansatz 

A(£) -> ar 1/2 , J(cj, I) -> £- 1/2 J{ujV1) 

which leads to the following non-local differential equa- 
tion for J(y) with the scale-invariant variable y = uj^/J: 

d y J(y) = -4J(y) ((y - a) - 2 f° dy'^1 _ * \ , 

v Jo v-v 4 ?/ y 

(55) 

where t = 1 — 4 J °° dyJ(y). The spectral function de- 
pends on the parameter a and the boundary condition is 
given by a = 4 J °° dyyj(y). 

With the similar ansatz for the observable flow, i.e. 

h{£) = br 1/4 - s ' /4 , S(uj,£) -> t- 1 ^''/ i S(wV£) , 
we obtain the following non-local differential equation for 



S(y) with y = u)\i: 

d y S(y) = -2S(y) ((y - a) - 2 f° dy'^1 *\ 

j o u y y 

+ 2J(y)(b + 2 r d y'^l + -) , (56) 

v Jo y-y y ' 

where t' = (s'+t)/2 and t" = J °° dyS(y). The parameter 
b is given by b = -4 / Q °° dyyS(y)/(s' + 1). 

The above ansatz assumes that a = A(^)^ 1 / 2 is con- 
stant. This is not the case, though - for the entire 
regime that is accessible by our numerical integration. 
We thus extend the asymptotic behavior of the tunnel- 
matrix element A(£) by setting A(£) = a{£)l~ 1/2 with 
a(£) — > a. The generalization induces a shift in the spec- 
tral function, i.e. J(uj,£) -> £~ 1/2 J(ujVI + a{£)) with 
[a — a(£)]/a(£) = const. 

The same extensions are also necessary in the case 
of the asymptotic observable flow, i.e. h(£) — > 
6 (^-i/4-«74 with b ^ b and S{uj,£) -> 

l-i/4-s' l±S{uj\T£ + a{£)) with [b-b{£)\/a{l) = const. All 
limits hold for t — > oo. A numerical fit for Ohmic dis- 
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FIG. 8: Left hand side: The asymptotic spectral functions 
S(y) obtained from the extended flow for Ohmic coupling with 
a = 0.1 and lj c /Ao = 10 at various values of £. Right hand 
side: The "analytic" solution following from Eq. (1581 for the 
corresponding values of b{£) = h^l 1 ^ 2 as they follow from the 
flow equations. A numerical fit suggests b — > 2.0 for £ — » oo. 



With At" = —s'b/a which follows from the boundary con- 
dition, the "analytic" solutions are again in good agree- 
ment with the flow equation solutions after a shift by 
the corresponding a(£). This is demonstrated in Fig. [SI 
where on the lhs the flow equations were integrated nu- 
merically and on the rhs Eq. Ij58(l was used. The various 
values of b(£) were obtained from the numerical integra- 
tion of the flow equations. 

Keeping in mind that J(y) is exponentially small for 
y — > 0, we are allowed to neglect the inhomogeneous 
contribution in Eq. I|56() or (|58|l in this regime. The 
asymptotic behavior of S(y) is then given by S 2 (y) oc 
J(y)y s ' for y - 0. With J(y) -» J and S(y) - 5 
for y — > 0, we thus obtain the £- independent correlation 
function for u) — * 



Jo 



(59) 



sipation with a = 0.1 and w c /A = 10 - based on the 
resulting differential equation for a{£) and b(£) and on 
the fact that [a — a(£)]/[b — b{£)\ — const - suggests that 
a = 4.5 and b = 2.0. 

But in the following, we want to demonstrate that 
the ^-dependent shift in the spectral function does not 
essentially alter the functional dependence of J(y) and 
S(y). For this, we define J(y) = J(y + a). Substi- 
tuting y — > y + a in Eq. (|55J) and setting a — * oo on 
the lower bound of the integral, the boundary condition 
yields t — and we obtain 

d y J{y) = -AyJ{y)(l-A H d y'^\). (57) 

v Jo y-y J 

This is the same non-local differential equation as in the 
case of the dissipative harmonic oscillator and the spin- 
boson model with form-invariant flow. But this time the 
parameter a of Eq. <|24|l is given by a = 2. The "analytic" 
result of the above equation is in good agreement with 
the flow equation result after a shift by the corresponding 
a(£). This is demonstrated in Fig. [7J where on the lhs 
the flow equations were integrated numerically and on 
the rhs Eq. H57J1 was used. The initial spectral function 
was taken to be J(lu) = 2aa->0(a-> c — lu) with a = 0.1 and 
uj c / Ao = 10 and the various values of a(£) were obtained 
from the numerical integration of the flow equations. 

The same transformation can be done for the asymp- 
totic function of the observable. With S(y) = S(y + a) 
this yields the following equation: 

d y S(y) = -2S(y)(y - 2 j°° dy> ^1 + jr^—) 
V J-oo y-y 4(y + a)/ 

+ 2J(y)(b + 2 f°° dy'^l + -^—) (58) 

v J~oo y-y y + aJ 



For an initial spectral function J (to) oc u> s , s — s' is 
verified by the numerical data and we thus obtain the 
correct low-energy behavior C{to) oc lu* . 

We want to note that the equations for the asymptotic 
behavior hold independent of the initial coupling strength 
and of the bath type. The initial conditions only enter 
via the parameters a and b for which the corresponding 
asymptotic functions are obtained. But we also want to 
mention that for higher coupling strength a and for lower 
values of s the asymptotic regime = X z a is only reached 
for rather large £. For an initial spectral function with 
a = 0.1, lo c /A = 10, K/A = 1 and s = 0.3 we have 
e.g. A^*>10 20 , £* denoting the asymptotic regime were 
= A^ holds on all energy scales. The corresponding 
parameters at I* are given by a = 2.23 and b — 1.43. 

Nevertheless, one can extend the range of the asymp- 
totic regime if one drops the assumption A^ = \ z a , i.e. 
one makes the ansatz J k (uj,£) — > £~ 1 l 2 J k {w\f£) and 
S k (w,£) -> £- 1 / i - s '/ i S k {uVI) for k = y,z . The result- 
ing set of non-local differential equations for the functions 
J k (y) recovers the flow equation results also for £*<£*, 
if the second boundary condition is suitably chosen. In 
fact, £* can be chosen such that £* <C £*■ 

To conclude the discussion, we now derive the asymp- 
totic conservation law, in analogy to Eq. I|31() . For this, 
we define the ^-dependent correlation function C(lo,£) = 
J2a Xa(V( 1J ^ u a) an d the asymptotic spectral function 
C a (cj,£) = C(lu) - C(lu,£). For £* being in the quasi- 
asymptotic regime, we obtain from the above flow equa- 
tions C a {uJ,£*) = lo s ' I(y*) where I{y*) is given by 
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I(y*) = i I dyy- s '[bSy(y) + S z (y) dy' S* (y')G + (y,y') + (y) dy' S z (y')G- (y,y') + t" — (60) 



We defined G ± (y,y') = (y - y')' 1 ± {y + y')- 1 and y* = 

uiVJ*. The above expression simplifies considerably for 
£* being in the asymptotic regime^* 

3. Numerical results 

We will now present the numerical results for the equi- 
librium correlation function C(lo). 

On the left hand side of Fig. the spectral func- 
tion C{u)/lo s is shown for various bath types s. The 
initial coupling function was chosen to be J(cj) = 
2aK 1 - s Lj s Q{uj c - uj) with a = 0.1, K/A = 1, and 
u! c / Aq = 10. As can be seen, C(u>)/ui s is constant for 
low energies. But notice that d^C \u>) / 'uj s ^ for ui — > 0. 
This point was already mentioned in the discussion of the 
dissipative harmonic oscillator (see Fig. Q}. 

Apart from the sum rule (a z (£)) « 1 for all £, another 
criterion for the assessment of the quality of the flow 
equations is given by the generalized Shiba relation^ 

lixu^ = 2«*W r^) 2 (61) 

The results are given in TableQ] So far, to our knowledge, 
the generalized Shiba relation for sub-Ohmic dissipation 
can only be checked by the above procedure. 

For Ohmic dissipation, the outlined approach can be 
contrasted with previous results. As initial spectral func- 
tion we use J(u)) = 2au}<d(uj c — u>) with a = 0.1 and 
w c / Aq — 10. On the right hand side of Fig. the spec- 
tral function C(uj)/uj is shown in comparison to the result 
obtained by integrating the flow equations of Sec. IIIA 
with f(w,£) = -(uj-A(£))/(u + A(£)). It is obtained by 
integrating the flow equations up to Aj^* = 100 and then 
using the conservation law of Eq. (|17fl . For higher values 
of A§£*>100, the spectral function shows instabilities for 
energies close to the resonance 10 w A*. The situation 
becomes worse for higher values of a and generally for 
sub-Ohmic baths. For a detailed discussion, see Ref. 0. 

In Tabic ITT1 the two flow equation approaches are con- 
trasted with the help of the generalized Shiba relation for 
various weak Ohmic couplings. Throughout the parame- 
ter regime, flow equations with universal asymptotic be- 
havior yield better results. For a = 0.2, the non-universal 
approach yields no stable result in order to check the gen- 
eralized Shiba relation. 

On the right hand side of Fig. C(ui)/uj is also shown 
as obtained from the non-interacting blib approximation 
(NIB A) for suitable effective tunneling A c ff. 19 The curve 
coincides with the one obtained from the universal ap- 
proach for short and intermediate time scales. For long 



r 





lim^o C(uj)/uj s 


2aK 1 - s (f™ duC(u)/uy 


Ratio 


s=1.0 


0.364 


0.391 


0.93 


s=0.8 


0.356 


0.390 


0.91 


s=0.5 


0.384 


0.432 


0.89 



TABLE I: The generalized Shiba relation for various bath 
types s with a = 0.1, K/A = 1, and u c /A = 10. 



R 


a = 0.01 


a = 0.05 


a = 0.1 


a = 0.2 


Universal 


1.00 


0.97 


0.93 


0.86 


Non-Uni. 


1.03 


1.19 


1.38 





TABLE II: Contrasting the universal and non-universal flow 
equations with the help of the generalized Shiba relation, 
R = [lim^oCHH/pa^ 00 dtuC{uj) /tuf], for a Ohmic 
bath with cj c /Ao = 10 and various coupling strengths a. 

time scales, the graph diverges as C(uo)/ui oc uj~ 2a , which 
is incorrect. 

The universal flow equation approach is applicable for 
not too small bath types s > s* and not too large 
coupling strengths a < a*. As criterion we again 
use the generalized Shiba relation and define the ratio 
R = [Hm^oCi^/ui^/fiaK 1 - 8 ^ dujC{oj)/uj) 2 }. We 
demand that |1 — i?|<0.3. For Ohmic coupling and 
u} c /Ao = 10, e.g., we have a*<0.4 and for a fixed cou- 
pling strength, a — 0.1 say, we have s*>0.3 (K/ Aq = 1). 
These parameter regimes are further discussed in Ref. l24l 
where also the case of finite bias is treated. 

If one is only interested in the effective tunnel-matrix 
element or the half-width of the resonance, the flow equa- 
tions with universal asymptotic behavior can yield the 
correct result even if the Shiba relation is not fulfilled 
satisfactorily. 



IV. CONCLUSIONS 

In this work, we applied the flow equation method to 
two dissipative systems, the exactly solvable dissipative 
harmonic oscillator and the spin-boson model. Investi- 
gating the former model, we have pointed that there are 
two basic universality classes for flow equations, i.e. ei- 
ther Woo > and Eq. (|10|) holds or w^, = 0. The for- 
mer class was labeled non-universal since the asymptotic 
behavior depended on the initial conditions of the sys- 
tem. The latter class was labeled universal since the fixed 
point Hamiltonian was given by the non-interacting bath, 
only. The universal asymptotic behavior of the spectral 
coupling function as well as of the observable are charac- 
terized by scale- invariant, non-local differential equations 
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s=1.0 




a s=0.8 


1 


I s=0.5 

I 







Universal 

Non-Uni. 

NIBA 




FIG. 9: Left hand side: The spectral function C(u>)/uj a for 
J(u>) = 2aK 1 ~ s Lu s e(u Jc - uj) with a = 0.1, K/Ao = 1, 
w c /Ao = 10 for various bath types s. Right hand side: The 
spectral function C(ui)/u> for J(uj) = 2aLuQ{uj c — to) with 
a = 0.1 and oj c /Ao = 10 as they follow from universal (solid 
line) and non-universal (dotted line) asymptotic behavior (for 
details, see text). The NIBA result is shown as dashed line. 



from which the low-energy behavior of correlation func- 
tions can be deduced analytically. 

These ideas were also applied to the non-trivial spin- 
boson model. It was mentioned that the form-invariant 
flow of the Hamiltonian necessarily displays universal 
asymptotic behavior when the reflection symmetry is bro- 
ken. The resulting flow equations, though, did not yield 
the correct normalization condition nor the correct low- 
energy behavior of correlation functions. Nevertheless, 



an effective energy scale was recovered which was in good 
agreement with the predictions of the NIBA. 

Due to the obvious shortcomings of the form-invariant 
flow, we set up extended flow equations which also dis- 
played universal asymptotic behavior. The resulting flow 
equations yielded very good results for a wide range of 
the parameter regime where also sub-Ohmic baths could 
be included in the treatment. We further deduced the 
correct low-energy behavior of correlation functions from 
the universal asymptotic behavior and showed that the 
normalization condition and the Shiba relations were sat- 
isfied within numerical errors for a certain parameter 
regime. 

Universal asymptotic behavior is crucial in order to 
determine correlation functions within the flow equation 
approach since only then one can assure that the domi- 
nant weight of the correlation functions is accounted for 
by the stable intermediate flow. Apart from this tech- 
nical preference, the trivial fixed point Hamiltonian also 
accounts for the correct spectrum which is not super- 
posed by the spectrum of the isolated system. We also 
want to mention that the asymptotic flow of the spectral 
function is described by the same non-local differential 
equation - independent of the model or the truncation 
scheme, but only differing by one parameter. 

The presented procedure can easily be extended to 
other dissipative models by starting from the canonical 
generator r\ c = [Hb , H ] where Hb denotes the non- 
interacting bath. 

It is a pleasure to thank A. Mielke for innumerable 
discussions on the flow equation method. This work was 
funded by the Deutsche Forschungsgesellschaft. 
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